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ABSTRACT 

Substellar objects, like planets and brown dwarfs orbiting stars, are by- 
products of the star formation process. The evolution of their host stars may 
have an enourmous impact on these small companions. Vice versa a planet 
might also influence stellar evolution as has recently been argued. 

Here we report the discovery of a 8—23 Jupiter-mass substellar object orbiting 
the hot subdwarf HD 149382 in 2.391 days at a distance of only about five solar 
radii. Obviously the companion must have survived engulfment in the red-giant 
envelope. Moreover, the substellar companion has triggered envelope ejection 
and enabled the sdB star to form. Hot subdwarf stars have been identified as 
the sources of the unexpected ultravoilet emission in elliptical galaxies, but the 
formation of these stars is not fully understood. Being the brightest star of its 
class, HD 149382 offers the best conditions to detect the substellar companion. 
Hence, undisclosed substellar companions offer a natural solution for the long- 
standing formation problem of apparently single hot subdwarf stars. Planets and 
brown dwarfs may therefore alter the evolution of old stellar populations and may 
also significantly affect the UV-emission of elliptical galaxies. 

Subject headings: binaries: spectroscopic — stars: individual(HD 149382) — planetary 
systems — stars: low-mass, brown dwarfs — stars: horizontal-branch — galaxies: 
evolution 
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Introduction 



A long-standing problem in extragalactic astronomy is the UV-excess observed in the 
spectra of elliptical galaxies. This phenomenon is caused by a n old popula tion of helium 



Hob or 



20091 ). The origin 



burning stars, known as hot subdwarfs or sdBs (see review by 
of the UV-excess can, hence, be traced back to that of the sdB stars themselves. The 
formation of such stars remains a mystery as it requires an extraordinarily high mass loss 
on the red-giant branch. Hot subdwarfs often reside in close binaries, formed by ejection of 
the envelope of their red-giant progenitors through interaction with the stellar companion. 
However, for half of the known sdBs no such companions could be found, requiring a yet 
unknown sdB formation channel. 

After finishing core-hydrogen-burning the progenitors of sdBs leave the main sequence 
and evolve to red giants before igniting helium and settling down on the extreme horizontal 
branch. Unlike normal stars, the sdB progenitors must have experienced a phase of 
extensive mass loss on the red giant branch, in order to explain the high temperatures and 
gravities observed at the surface of hot subdwarf stars. After consumption of the helium 
fuel they evolve directly to white dwarfs avoiding a second red-giant phase. What causes 
this extensive mass loss remains an open question. 

The riddle of sdB formation is closely related to other long-standing problems regarding 
old stellar populations, which have been discussed for decades. The morphology of the 
horizontal branch in globular clusters, especially the existence and sha pe of its extre me hot 



part (Extreme Horizontal Branch, EHB), is still far from understood (jCatelan 



2009|). Hot 



subdwarfs are also regarded as the dominant source of the UV-excess in early type galaxies, 
where no active star formation is going on and hence no UV-emission from young massive 
stars is expected . Hot subdwarf formation is the key to understa nding the physics behind 



this phenomenon and a debate is going on whether single star (JYi 



20081 ) or binary evolution 
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( IHan et al.l 120071 ) explain the observed UV-excesses. 



About ha 



(IMaxted et al. 



of the sdB stars reside in clo se binaries with periods as short as ~ 0.1 d 



2001 



Napiwotzki et al. 



2004aj). Because the components' separation in these 



systems is much less than the size of the subdwarf progenitor in its red- giant phase, thes e 



2002 . 



syste ms must have experienced a common-envelope and spiral-in phase ( lHan et al. 
20031 ) . In such a scenario, two main-sequence stars of different masses evolve in a binary 
system. The more massive one will first reach the red-giant phase and at some point fill 
its Roche lobe, where mass is transferred from the giant to the companion star. When 
mass transfer is unstable, the envelope of the giant will engulf the companion star and form 
a common envelope. Due to friction with the envelope, the two stellar cores lose orbital 
energy and spiral towards each other until enough orbital energy has been deposited within 
the envelope to eject it. The end product is a much closer system containing the core of 
the giant, which then may become an sdB star, and a main-sequence companion. This 
companion evolves to a white dwarf after another phase of unstable mass transfer. 

The common-envelope ejection channel provides a reasonable explanation for the extra 
mass loss required to form sdB stars. But for about half of all analysed subdwarfs there is 
no evidence for close stellar companions as no radial velocity variations are found. Although 
in some cases main sequence companions are visible in the spectra, it remains unclear 
whether these stars are close enough to have interacted with the sdB progenitors. Among 
other formation scenarios, the merger of two helium white dwa rfs has often been suggested 



to explain the origin of single sdB stars (IHan et al 



20oJ . 



20031 ) . Merging should result in 



rapidly spinning stars, which is not consistent with observations. A recent analysis of single 
sdB stars revealed that their v ro t sin ^-distribution is consistent with a uniform rotational 



velocity v ro t ~ 8 km s 1 and randomly oriented polar axes (jGeier et al. 



2009|). 



The planet discovered to orbit the sdB pulsator V 931 Peg with a period of 1 170 d 
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and a separation of 1. 7 AU was the first 



phase of its host star (jSilvotti et al. 



planet found to have survived the red-giant 



20071 ). Serendipitous discoveries of two substellar 



companions around the eclipsing sdB bin ary HWVir at d istances of 3.6 AU and 5.3 AU 



with orbital periods of 3 321 d and 5 767 d (ILee et al. 



20091 ) and one brown dwarf around the 



similar system H S 0705+6700 with a period of 2 610d and a separation of less than 3.6 AU 



( jOian et al. 



20091 ) followed recently. These substellar companions to hot subdwarfs have 
rather wide orbits, were not engulfed by the red giant progenitor and therefore could not 
have influenced the evolution of their host stars. 



2. Observations and Analysis 

HD 149382 is the brightest core helium-burning subdwarf known. The first hint that 
this star could show very small radial velocity (RV) variations was fo und during our survey 



aimed at finding sdBs in long period binaries (jEdelmann et al. 



20051 ). We obtained 15 high 



resolution spectra (R = 30 000 — 48 000) within four years with three different high resolution 
spectrographs (ESO-2.2m/FEROS, McDonald-2.7m/Coude, CAHA-2.2m/FOCES). One 
additional spectrum obtained with ESO-VLT/UVES at highest resolution (R « 80 000) was 
taken from the ESO-archive. 

In order to measure the RVs with highest possible accuracy we fitted a set of 
mathematical functions (polynomial, Gaussian & Lorentzian) to all suitab le spectral lines 



with w avelengths from about 4 000 A to 6 700 A using the FITSB2 routine ( Napiwotzki et al. 



2004bj). The formal deviation along the whole wavelength range was 0.2 km s _1 at best. 
In order to check the accuracy of this measurements, we also obtained RVs from telluric 
and night sky lines and reached similar accuracies. Since telluric and night sky lines have 
zero RV we used them to correct the measured RVs for calibration errors. The applied 
corrections were usually below 1.0 km s _1 . Since we used four entirely different instruments 
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and obtained consistent results other systematic effects should be neligible. The RV 
measurements are given in Table [TJ 

The period search was carried out by means of a periodogram based on the Singular 
Value Decomposition (SVD) method. A sine-shaped RV curve was fitted to the observations 
for a multitude of phases, which were calculated as a function of period. The difference 
between the observed radial velocities and the best fitting theoretical RV curve for each 
phase set was evaluated in terms of the logarithm of the sum of the squared residuals (x 2 ) as 
a function of period. This method finally results in the power sp ectrum of the data s et which 



allows to determine the most probable period P of variability (lLorenz et al.l 119981 ). The 
formal significance of the best orbital solution (P = 2.391 ± 0.002 d, K = 2.3 ±0.1 km s -1 , 
7 = 25.3 ± 0.06km s _1 ) exceeds the 3a-limit (see Figured]) and a very small mass function 
f(M) = 3.8 x 10~ 6 results. The radial velocity curve is shown in Fig. [2j The formal 
probability that the next best alias periods at about 4.8 d and 6.0 d are correct is less 
than 5 %. Even if one of these longer periods should be the correct one, the mass function 
increases only by a factor of 1.5 at most and our conclusions still remain valid. 

The atmospheric parameters effective temperature T e g, surface gravity log g and 
helium abundance were determined by fitting simultaneously 17 hydrogen and helium 
lines in high resolution, high- S/N FEROS and UVES spectra with NLTE model spectra 



(the method is described in 



Lisker et al 



2005 



Geier et al 



20071 ). The parameters 



(T P ff = 35 5 ± 5 00 K, log g = 5.80 ± 0.05) are in good agreement with the result of 



Saffer et al 



(|l99J): T eff = 34 200 ± 1500 K, log g = 5.89 ±0.15. 



The mass of the unseen companion can be derived by solving the binary mass function 
f m = Mj? Q sin 3 z/(M comp + M sdB ) 2 = PK 3 /2nG. In order to obtain a unique solution, the 
mass of the sdB primary as well as the inclination of the system must be known. Due to 
the excellent quality of the data available for HD 149382, constraints can be put on both 
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crucial parameters. 



The distance to this star can be deri ved directly using a trigonometric parallax 



obtained with the HIPPARCOS satellite (Ivan Leeuwen 



20071 ). We derive the angular 



diameter by comparing the surface flux in the V band computed from a model atm osphere 



with the derived atmospheric parameters to the observed value ( iMermilliod 



19911 ) ■ Using 



the trigonometric distance we can deriv e the stellar radius, and from the surface gravity the 



mass of the sdB (IRamspeck et al 



200 lh 



Taking the uncertainties of all parameters into account (V magnitude, T e g, logg, 
parallax) the possible mass range for the sdB is 0.29 — 0.53 M Q . This range is consistent 
with the canonical ma ss of 0.47 M & derived from single and binary evolution calculations 



( iHan et al. 



2002 



20031 ). Without further constraints on the inclination only a lower limit 



to the mass of the unseen companion can be calculated. 

The minimum companion mass lies between 0.006 M and 0.01 M^, well below the 



stella r limit of 0.075 M Q to 0.083 M & depending on the metallicity (IChabrier fc Baraffe 
19971 ). which is the lower limit where core hydrogen burning can be ignited and a star can 
be formed. The lower the inclination of the binary is, the higher is the mass of the unseen 
companion. Assuming a random distribution of orbital plane inclinations , the probability 



19921 ) . The probability 



for a binary to fall below a certain inclination can be derived (iGrayl 
for the companion to have a mass of more than 0.08 M is just 0.8 %. The probability that 
the mass of the unseen companion exceeds the planetary limit of 0.012 M & defined by the 
IAuU is only 33 %. 



However, we can constrain the mass of HD 149382 b even further. Due to the very 



*Position statement on the definition of a "planet". Working group on extrasolar planets 
of the International Astronomical Union, http://www.dtm.ciw.edu/boss/definition.html 
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high resolution of the UVES spectrum the broadening caused by the projected rotational 
velocity of the star could be measured from the metal lines although it turned out to 
be as small as t> rot sini = 4.9 ± 1.4 km s _1 . In order to derive v TOt sini, we compared 
the observed high resolution (R = 80 000) UVES spectrum with rotationally bro adened 



synth etic line profiles. The profiles were computed using the LINFOR program (iLemke 
19971 ). A simultaneous fit of elemental abundance, projected rotational velocity and radial 
velocity was then perform ed separately for every ident ified line using the FITSB2 routine 



(INapiwotzki et al. 



2004bl ). The method is described in 



Geieretal 



(120081) . 



Assumi ng that HP 14938 2 rotates with the standard velocity of 8 km s 1 infered for 



single sdBs (jGeier et al. 



20091 ) the inclination can be constrained to i = 26° — 52°. The 



companion mass range is derived to be M2 = 0.008 — 0.022 M© = 8 — 23 Mj consistent 



with a gas giant planet or a low mass brown dwarf. Adoptin g the statistical 



2002 



y mos t likely 



20031 ) the 



inclination i = 52° and the canonical sdB mass of 0.47 M & (lHan et al. 
companion mass is 0.011 M & = 12 Mj, which places HD 149382 b just below the planetary 
limit. The separation of star and companion is 5 — 6R Q . All relevant measurements and 
parameters of the HD 149382 system are summarized in Table [2j 



3. Discussion 



When the progenitor of HD 149382 evolved through the red-giant phase, it expanded 
to a radius of ten times the present orbital separation. The initial separation must have 
been larger (about 1 AU) and the companion spiralled- in due to interaction with the giant's 
envelope until the envelope was ejected. Despite the very high lo cal tempera ture inside the 



envelope (300 000 to 400 000 K at 5 — 6 R Q from the giant's centre, 



Sokcr 



19981 ) the substellar 



companion survived. The compa nion of the sdB star AA Dor has also been suggested to be 



a brown dwarf in a 0.26 d orbit (IRauch 



2000 : 



Rucinski 



20091 ). This conclusion, however, 
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is rendered uncertain as IVuckovic et al.l (120081 ) derive a higher mass indicating that the 
companion is a star. 



Sokerl (119981 ) suggested that sub-stellar objects like brown dwarfs and planets may 
also be swallowed by their host star and that common-envelope ejection could form hot 
subdwarfs. Substellar objects with masses higher than rj 10 Mj were predicted to survive 
the common-envelope phase and end up in a close orbit around the stellar remnant, while 
planets with lower masses would entirely evaporate. The stellar remnant is predicted to 
lose most of its envelope and evolve towards the EHB. The orbital period and mass w e 



derive d for HD 149382 b is in excellent agreement with the predictions made by 



Soker 



(119981 ) . A similar scenario ha s been proposed to expla in the formation of apparently 



single low mass white dwarfs (iNelemans fc Tauris 



1998). The discovery of a brown dwarf 



with a mass of 0.053 ± 0.006 M in a 0.08 d orbit around such a white dwarf supports 
this scena rio and shows that substellar companions can influence the outcome of stellar 



evolution (IMaxted et al. 



20061 ). 



The discovery of planets and brown dwarfs around sdBs and especially the close-in 
substellar companion of HD 149382 may thus have important implications for the still 
open question of sdB formation. The extraordinary quality of the photometric data 



was a prerequisite f or the detection of t he subste 



and HS 0705+6700 (jSilvotti et al 



2007 



Lee et al- 



ar companions in V 93 1 Peg, HWVir 



2009 



Qian et al 



20091 ). Finding such 



companions orbiting three of the best observed sdBs cannot be mere conincidence and leads 
to the conclusion that substellar objects may often be associated with sdBs. HD 149382 
is the brightest sdB known. Hence the quality of the spectroscopic data is also very high. 
It is not easy to detect such small RV variations even in high resolution spectra. The fact 
that we found them in the case of HD 149382 leads to the conclusion that close-in planets 
or brown dwarfs may be common around apparently single sdB stars. They were just not 
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detected up to now. Hence all apparently single sdBs may have or had close brown dwarf 
or planetary companions, although those of lowest mass may have evaporated. 

HD 149382 b provides evidence that substellar companions can decisively change the 
evolution of stars, as they trigger extensive mass loss. They could be responsible for the 
formation of the single hot subdwarf population. These stars are not only numerous in our 
Galaxy, but also make elliptical galaxies shine in ultraviolet light. 

This letter was based on observations at the La Silla Observatory of the European 
Southern Observatory for programme number 077.D-0515(A). Some of the data used in 
this work were obtained at the McDonald Observatory. Based on observations collected at 
the Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max-Planck Institut fur Astronomie and the Instituto de Astrofisica de Andalucia (CSIC). 
Some of the data used in this work were downloaded from the ESO archive. S.G. gratefully 
acknowledges financial support from the Deutsche Forschungsgemeinschaft through grant 
He 1356/40-4. We thank Sebastian Miiller for reducing the archival UVES spectrum. 
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Fig. 1. — Power spectrum of HD 149382. — logx 2 is plotted against the orbital period in 
days. The region around the most prominent period is plotted in the small window. The 
formal significance exceeds the 3<r-limit. 
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Fig. 2. — Radial velocity curve of HD 149382. The plot shows the radial velocity plotted 
against orbital phase (diamonds: McDonald-2.7m/Coude, rectangles: CAHA-2.2m/FOCES, 
upside down triangles: ESO-2.2m/FEROS, triangle: ESO-VLT/UVES). The RV data was 
folded with the most likely orbital period. The residuals are plotted below. 
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Table 1: Radial velocities of HD 149382 



mid-HJD RV [km s *] Instrument 



2452497. 


.49150 


27 


.0 


± 





.2 


FEROS 


2452497. 


.51327 


26. 


.9 


± 


0. 


.2 


FEROS 


2452497. 


.56524 


26. 


.9 


± 





.2 


FEROS 


2452891. 


.30690 


25. 


.3 


± 





.2 


FOCES 


2452892. 


.30170 


25. 


.7 


± 





.9 


FOCES 


2452893. 


.32160 


23 


.9 


± 





.6 


FOCES 


2453784. 


.83294 


23. 


.8 


± 





.2 


UVES 


2453904. 


.73904 


22. 


.9 


± 





.2 


FEROS 


2453931 


.76614 


26. 


.8 


± 





.2 


Coude 


2453932 


.71806 


26. 


.1 


± 





.2 


Coude 


2453932 


.72485 


25. 


.7 


± 


0. 


.2 


Coude 


2453932 


.74337 


25. 


.9 


± 


0. 


.2 


Coude 


2453932 


.74882 


25. 


.7 


± 





.2 


Coude 


2453932 


.84790 


24. 


.8 


± 





.2 


Coude 


2453932 


.85560 


25. 


1 


± 





.2 


Coude 


2453986. 


.54830 


24. 


.8 


± 





.2 


FEROS 
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Table 2: Parameters of the HD 149382 system 



Trigonometric parallax a 


7T 


[mas] 


13.53 ±1.15 


Distance 


d 


[P C ] 




Visual magnitude 13 


my 


[mag] 


8.947 ±0.009 


Atmospheric parameters 


Effective temperature 


T e s 


[K] 


35 500 ±500 


Surface gravity 


log# 




5.80 ±0.05 


Helium abundance 


logy 




-1.44 ±0.01 


Projected rotational velocity 


v TO t sin % 


[km s _1 ] 


4.9 ± 1.4 


Orbital parameters 


Orbital period 


p 


[d] 


2.391 ±0.002 


RV semi-amplitude 


K 


[km s _1 ] 


2.3 ±0.1 


System velocity 


7 


[km s _1 ] 


25.3 ±0.06 


Binary mass function 


f(M) 


[M ] 


3.8 x 10~ 6 


Derived parameters 


Sub dwarf mass 


M sdB 


[M @ ] 


0.29-0.53 


Orbital inclination 


i 


[°] 


26 - 52 


Companion mass 


■^comp 


[Mj] 


8-23 






[M & ] 


0.008 - 0.022 


Separation 


a 


[Re] 


5.0-6.1 



a The trigonometric parallax was taken from the new reduction of the HIPPARCOS data (jvan Leeuwenl l2007). 
h The visual magnitude is taken from iMermilliodl (|1991l ). 



